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Abstract— : In this research, we formulated the equation for water wave 
velocity potential on sloping bottoms. The resulting velocity potential 
equation for sloping bottoms closely mirrors that of flat bottoms, 
simplifying its application in sloping terrain scenarios. By exploring 
velocity potential on sloping bottoms, we derived conservation equations 
governing wave constant changes as waves transition from deep to 
shallow waters. These equations encompass the wave number 
conservation and energy conservation principles. Utilizing these 
conservation equations, we developed a comprehensive wave 
transformation model. The one-dimensional model focused on shoaling 
and breaking phenomena, while the two-dimensional model delved into 
refraction-diffraction, shoaling, and breaking. The model's framework 
allows for straightforward extensions, facilitating future advancements in 


the field. 


I. INTRODUCTION 

Waves undergo significant transformations as they travel 
from deep waters to shallow coastal areas. These changes 
involve alterations in their fundamental properties, such as 
a decrease in wavelength and an increase in wave 
amplitude, eventually leading to the breaking of waves. 
Consequently, a _ velocity potential capable of 
accommodating these variations in wave constants 
becomes essential. 
The velocity potential of water waves is typically 
determined by solving the Laplace equation. Dean (1994) 
utilized a flat bottom approach to formulate this potential, 
resulting in constant wave number and wave amplitude. 
However, this method obscures the detailed changes in 
wave constants. To address this limitation, researchers 
have explored the wave group concept, focusing on the 
conservation of energy within wave groups to model 
variations in wave number and wave height due to changes 
in water depth (Dean, 1994). Several wave transformation 
models have emerged from this concept, with the Mild 
Slope Equation (MSE) pioneered by Berkhoff (1972) 
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being the most renowned. Subsequent researchers such as 
Booij (1983), Davies and Heathershaw (1984), and Porter 
(1995) have further developed and refined the MSE. Hsu, 
Lin, and Ou (2006) introduced the Complementary Mild 
Slope Equation (CMSE), while Lan, Yuan-Jyh, Hsu, Tai- 
Wen, Lin, Ta-Yuan, and Liang, Shin-Jye (2012) proposed 
a wave transformation model utilizing the Higher Order 
Mild-Slope Equation. 

In our research, we delved into solving the Laplace 
equation within the context of a sloping bottom. This 
investigation enabled us to uncover the distinctive 
alterations in wave constants, specifically the conservation 
laws embedded within the solution of the Laplace equation 
at the sloping bottom. These fundamental conservation 
laws serve as the foundation upon which various wave 
transformation models can be constructed, such as shoaling 
breaking and diffraction refraction. 


Il. GENERAL SOLUTION OF LAPLACE’S 
EQUATION 


The Laplace equation is, 
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ap , ap _ 
dx2 az2 = 0 ..(1) 


o(x,z,t) represents the potential velocity, where x 
represents the horizontal axis, z represents the vertical axis, 
and t represents time. The solution of this equation 
employs the separation variable method, as outlined by 
Dean (1994), which involves the application of periodic 
lateral boundary conditions. In the separation variable 
method, the velocity potential @(x, z, t) is conceptualized 
as the product of three distinct functions, namely, 


p(x, z,t) = XOTE) 

X(x) is a function of x only, Z(z)is a function of z only, 
and T(t) is a function of t only. The properties of these 
functions are crucial to remember when formulating 
various equations using the velocity potential equation 
obtained through the separation variable method of the 
Laplace equation solution. 


Dean (1994) obtained the general solution of the Laplace 
equation in the form of, 


o(x,z,t) = (A cos kx + B sin kx ) 
(Ce? + De™)sinot ...(2) 


Where, 

X(x) = (Acos kx + B sinkx ) 

Z(z) = (Ce*? + De™*7) 

T(t) = sin ot 

k and ø are wave constants, where k represents the wave 


number k = = while L being the wavelength and o 


represents the angular frequency o = = and T being the 
wave period. 
On the other hand, A, B, C, and D are constants of the 
solution that still need to be determined. To derive the 
values of these constants, the equation is manipulated at 
characteristic spatial points, where cos kx = sin kx, 
consequently, (2) can be expressed as, 
(x,z,t) = (A+B) coskx (Ce + De-*“) sin ot 
The solution constants will be formulated by applying the 
bottom kinematic boundary condition as follows. 
dh 
Yoh Tx 


Wp bottom vertical water particle velocity at z = —h 


W_ph= 


u_,, bottom horizontal water particle velocity at z = —h 


h water depth towards still water level. 
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dh ; ; 
T Slope of the water, negative for waves propagating 


from deep water to shallow water.. 


Vertical water particle velocity is, 


d 
w(x,z,t) = Ti 


= —(A +B )k coskx (Ce — De~*) sin ot 


w_, = —(A +B )k cos kx (Ce™™™" — De*") sin ot 
Horizontal water particle velocity is, 


d 
u(x, z,t) = eL 


= (A+B )k sinkx (Ce*? + De~”) sinot 


Uu_n = (A+B )k sinkx (Ce™*" + De*") sinot 


For (A+B) + 0 and at the characteristic point where 
cos kx = sinkx and sinot #0, the bottom kinematic 
boundary condition equation becomes, 


(Ce7kh — Dekh) = (Ce-** + Dekh) < pii (3) 


2.1 At the Flat Bottom 
At the flat bottom, where zZ = 0, (3) transforms to 
(Ce-** — Dek") = 0 
Obtaining 
C = De?” 
Being substituted to (2) 
@(x,z,t) = (A + B)coskx(De?" e" + De") 


sin(ot) 


@(x,z,t) = 2(A + B)De*"coskx 
ek(htz) okz 4 pk(h-z) 
2 
@(x,z,t) = 2(A + B)De*"coskx cosh k(h + z) sin (ot) 
Defined as A = 2A and B = 2B, thus 
@(x,z,t) = (A+ B)De""coskx cosh k(h + z) 


sin(ot) 


sin(ot) 
Is defined by the wave constant, 
G = (A+ B)De*r 
@(x, z,t) = G coskx coshk(h + z) sin(ot) (4) 


It is important to note that in G, there is an addition of two 
constants, namely (A + B), where Hutahaean (2003a) 
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demonstrated that A = B, thus G in equation (4) has a 
double value. 


2.2 At the Sloping Bottom 
Kinematic bottom boundary condition at the sloping 
bottom is expressed as, 


dh 
(Ce~*" — De”) = (Ce ™™ + Des) Tx 


dh 
C = De? Tay 
,— ah 
dx 
ah 
Defined as: a = —% 


C=Dae0e™ 

@(x,z,t) = (A+ B)coskx(D a e?*" ek + De") 
sin(ot) 

@(x,z,t) = 2(A + B)De*"coskx 


(< ek(n+z) +4 e~k(h+z) 


5 sin(ot) 


P(x, z,t) = 2(A + B)De*"coskx B(z)sin(ot) 
a ek(h+z) 4 e—k(h+z) 


Where : B(z) = 


2 
By defining the new constant G as stated in the previous 
section, the equation for the velocity potential on a 
sloping bottom is obtained as follows, 


@(x,z,t)=Gcoskx B(z)sin(ot) aaa (5) 


Subsequently, a comparative analysis was conducted 
between the values of p(z) and coshk(h +z), for 


f dh 
different bottom slope values represented by T when 


k(h+z) =O where 0 = 2.1, and tanh@8r = 1, as 
indicated in Table (1). The use of tanh Orm = 1 is specific 
to deep water conditions, a concept that will be elaborated 
upon in section (3). Consequently, 0 will be denoted as the 
deep water coefficient henceforth. 

aePT4e-9T 


Table (1). The comparison between B (0T) = z 


and cosh 0r 


dh s 
dx a B(On) | cosh 8r (%) 


-0.01 | 0.98 | 359.315 | 366.573 1.98 
-0.02 | 0.961 | 352.198 | 366.573 | 3.922 
-0.03 | 0.942 | 345.22 | 366.573 | 5.825 
-0.04 | 0.923 | 338.376 | 366.573 | 7.692 
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-0.05 | 0.905 | 331.662 | 366.573 | 9.524 
-0.06 | 0.887 | 325.075 | 366.573 | 11.321 
-0.07 | 0.869 | 318.611 | 366.573 | 13.084 
-0.08 | 0.852 | 312.266 | 366.573 | 14.815 
-0.09 | 0.835 | 306.039 | 366.573 | 16.514 
-0.1 | 0.818 | 299.924 | 366.573 | 18.182 


B(@m)-cosh On 
cosh 61 


Note : € = x100% 


In Table (1), the disparity between the values of B(@7) = 


aeFM4e—On ba 
and cosh @z, where tanh 8r = 1., is illustrated. 


hoon ath . : : dh 
It is evident that the difference increases with larger [Alor 


steeper bed slope. However, it is common for the seafloor 
to have a gentle slope with a bed slope of less than 0.01. 
Consequently, it can be inferred that f (@7) ~ cosh 67) or 
in other words, p (O71) can be approximated by ccosh Oz. 


Bi (Or. 


Table (2) Presents the comparison between BOn and 
tanh On at On = 2.1. 
dh B, (O71) 
dx a B(@n) | tanh 6x 
-0.01 0.98 1 1 
-0.02 0.961 1 1 
-0.03 0.942 1 1 
-0.04 0.923 1 1 
-0.05 0.905 1 1 
-0.06 0.887 1 1 
-0.07 0.869 1 1 
-0.08 0.852 1 1 
-0.09 0.835 1 1 
-0.1 0.818 1 1 
Subsequently, it is defined that 
a ek(n+z) _ e~k(ntz) 
B,(Z) = 5 
This equation is similar to sinhk(h+z), and a 
comparison is made between Ca and tanh @z for the 


deep water coefficient 0 = 2.1, as shown in Table (2). 


Table (2) illustrates that at 0 = 2.1, it is observed that 
Ba (67) 
B(6n) 
concluded that the coshk(h+z) function can be 
employed as a substitute for p(z). Consequently, the 


= tanh 6z. Based on this comparison, it can be 


velocity potential equation on a sloping bottom can be 
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replaced with the velocity potential equation on a flat 
bottom, namely equation (4), 


o(x,z,t) = Gcoskxcoshk(h+z)sinot ...... (4) 


Where the complete equation is, 
G 
o(x,z,t) = 3 (cos kx + sin kx) coshk(h + z) sin øt 


.(5) 
II. EQUATIONS OF CONSERVATION. 
3.1 Conservation Equation of Wave Number. 

In solving the Laplace equation using the separation 
variable method, it is assumed that the velocity potential is 
a product of three functions. In equation (4), where Z(z) 
represents a function of z only, it is defined as, 


Z(z) = coshk(h + z) 
Therefore, on a sloping bottom where changes in water 
depth h and wave number k occur, the following applies, 


daZ(z) 
=0 

dx 
sah RA dene) _ 5 
sinh k( 5 ae aa 
Where 
dk(h +z) 
=o 

dx 


In this equation, it is essential to determine the appropriate 
value for z . Hutahaean (2022a, b), in formulating the wave 
amplitude function using the kinematic free surface 
boundary condition, determined that the suitable value for 
; A f . 
ZiSZ=>. Therefore, the conservation equation of wave 
number becomes, 
A 

dk(h+5) 


dx 


=O 0 (6) 


Consequently, the equation obtained is, 
A 
k (n + 5) = constant 


These constants apply throughout the entire water body, 
including deep waters where 
A 
tanhk(h+5) =1 
For 


A 
k(h+5) = 6n 
The parameter 0 will be henceforth referred to as the deep 
water coefficient. Hutahaean (2023), through the study of 
the breaker depth index H, represents the breaker 
height and hp represents the breaker depth, recommends 
the value 0 = 2.1. 


Where small amplitude is assumed, (6) can be expressed is 
dkh 


—=0 
dx 
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For a wave moving from point x with water depth h, to 
point x + dx with water depth h,+5,, where hy+5, is 
shallower than hy, there exists a relationship. 

kyhy 


hy+5x 
It can be observed that k,45, is greater than k,, or the 


kx+6x = 


wavelength Ly+sx is smaller than Ly. This indicates a 
shortening of the wavelength when waves move from 
deeper water to shallower water. This means that in the 
wave number conservation equation, there is a 
phenomenon of shoaling, signifying a change in 
wavelength. 


As a consequence of equation (6), for small amplitudes, the 
following applies, 
d tanh kh 1 dkh 
dx ~ cosh? kh dx — j 
Or 
tanh kh = constant 


As known, the wave constants k, G and A are functions of 
water depth h, meaning they change with depth or are 
functions of h. Therefore, waves moving from deep water 
to shallow water will experience changes in these three 
wave constants. The water depth h is a function of the 
horizontal axis x, denoted as h = h(x). Consequently, k = 
k(x), G = G(x) and A = A(x). This implies the existence 


dk dG dA : . 
of values for TET and To The governing equation 
X X x 


controlling the changes in these wave constants is the 
continuity equation, where, with the variable velocity 
potential @, this equation becomes the Laplace equation. 
In deep water, tanh kh = 1. Therefore, for all water 
depths, including shallow water, tanh kh = 1. holds true. 
This condition is different from linear wave theory, where 
tanh kh = 1 changes with varying water depth. 
3.2 Conservation of Energy Equation 

As it is known, the wave constants k, G and A change 
with depth or are functions of water depth h. Therefore, 
waves moving from deep water to shallow water will 
experience changes in these three wave constants. The 
water depth h is a function of the horizontal axis x, denoted 
as h = h(x). Consequently, k = k(x), G = G(x) and A = 


: : dk dG dA 
A(x) leading to the existence of values for Te ae dan — 


The governing equation controlling the changes in these 
wave constants is the continuity equation, where, with the 
variable velocity potential @ this equation becomes the 
Laplace equation. 

Equation (5) is substituted into (1), and when evaluated at 
the characteristic point, it yields, 

a?G 


Ss =0 eee eee (7) 
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If the velocity potential given in equation (4) is used, and 
considering equation (7), 


GEER SG  — isons (8) 
dx dx 


Considering that G represents the rate of energy transfer 
per unit time, equation (8) can be referred to as the energy 
conservation equation. 


IV. THE WAVE AMPLITUDE FUNCTION 


By integrating the kinematic free surface boundary 
condition with respect to time, using the velocity potential 
given in equation (4), and evaluating it at characteristic 
points in space and time, Hutahaean (2022a,b) derived the 
water surface elevation equation as follows: 


n(x, t) = A cos kx cos ot 
A is wave amplitude. 


Gk 
20/2 


(1 — 5 cosh@m > n (9) 


This equation is referred to as the wave amplitude function. 
Yz is the coefficient in the time derivative of the Taylor 
series, due to the truncation of the series, to the first 
derivative only. Hutahaean (2022a) obtained y, = 1.4. 0 
is the coefficient for deep water. 


Taylor series for function f = f(x, t) 


df df 
f(x + 6x,t + ôt) = f (x,t) + bts + ôx 


dx 
ôt? a? a? ôx? a? 
aS a etoxi A q TT 
2! at? dtdx 2! dx2 


When truncation is performed, only the first derivative is 
used, Hutahaean (2022a), 


f(x + dx,t + dt) = f (x,t) + dt + ôx 
While for function f = f(x, z,t) 
df 


f(x + 6x,x + dx,t + dt) = f(x,z,t) + y3ôt =r 


+x L + 4622 
dx dz 


Equation (9) can be written as the equation for G as 
follows. 


20V2A 


ae k(1—*4) cosh 6m 


V. SHOALING-BREAKING MODEL 
In this section, a one-dimensional wave transformation 
model is developed, including only shoaling and breaking. 


Waves transitioning from deep water to shallower depths 
undergo alterations in the parameters G, k and A - a 
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phenomenon known as shoaling. During shoaling, the 
wave number k increases, leading to a shorter wavelength, 
and the wave amplitude A rises until reaching a breaking 
point at a specific depth. Developing shoaling and breaking 
models becomes straightforward by employing two 
conservation equations and the wave amplitude function. 


In the energy conservation equation, alterations in G 
interact with changes in k, while the wave number 
conservation equation shows a connection between 
variations in wave number k and changes in wave 
amplitude A due to alterations in depth h. Consequently, 
working with these conservation laws reveals an intricate 
interplay between the three wave constants: G, k and A 
Equation (8) can be written as, 

1dG_ = 1 dk 

Gdx  2kdx 

This equation is multiplied by dx and integrated, yielding: 


~ dG _ L peto dk 
x G 2x k 


1 
In Gx 45x —In Gy = 2 (n Retox —In ky) 
Or 


1 
en Gy-Z(n Ky45x—In Kx) 


wee chs 


This equation represents the relationship between the 


Gx+ôx = 


change in G concerning the change in k, for waves moving 
from point k, to point x + 6x where hy+5, < hy. From this 
equation, it is evident that knowing x + dx where hy45, < 
h, is necessary to calculate G,,5,. Therefore, what is 
needed is the equation describing the change in wave 
number k first. 


Equation (9) is differentiated with respect to the horizontal 


axis x 
dG _ 20, dA 
dx k (1 — KA) cosh 61 dx 
20Y2A dk 
k2 (1 — KA) cosh 61 dx 
a(1—"4) 
In the differential, the term is zero ~~ based on the 


wave number conservation equation. Considering (9), this 
equation can be written as, 

dG GdA Gdk 

dx Adx kdx 

Upon substitution into (8), 


eee —CUNsanteement (12) 


2kdx  Adx 


The wave number conservation equation can be written as, 
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dA o 2 m ‘) dk dh 

dx a 2/dx “dx 

Substituting into (12), 

dk 2k dh 

Peli (an 4) ax sagan (13) 


This equation is the change in wave number k equation that 
satisfies the energy conservation equation, wave number 
conservation equation, and kinematic free surface 
boundary condition, where (9) is derived from the 
kinematic free surface boundary condition equation 


In a very small interval, small 5x k,,, 5, can be assessed as, 
k ky + 6. = 
= x — 
x+6X x dx 
Furthermore, the change in wave amplitude A can be 
obtained by differentiating the wave amplitude function 
with respect to the horizontal axis x, 
dA G dk 
dx Aoy2 ax 


(1 = 5 cosh(On) n 14) 


where Z is obtained from (13). In a very small interval ôx, 


Ax+5x can be calculated as follows, 


dA 
Ax+ôx = A, + 6x dx 


In the following section, an example of shoaling model 
results is presented for waves with a wave period T = 7.15 
seconds and wave amplitude in deep water \Ay = 1.2 m 
meters, in waters with a bottom slope = = —0.01. The 
calculation parameters used are the deep water coefficient 


zZ = —0.01. The model results are presented in Figure 1. 


The wave period is calculated based on the input wave 
amplitude (Hutahaean (2023a)). 


2 Y3 i 
T= salar) 4 (sec) n (15) 
where y3 = 1.8 (Hutahaean (2022a)). 
12 
10 
€ 8 
6 
ð 
ka a 4 
2 
0 
0 2 4 6 8 10 
h (m) 


wave height H wavelength L 


Fig (1) The outcomes of shoaling-breaking model. 
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The results of the shoaling and breaking model (Fig (1)) 
show that breaking occurs at a breaker depth h, = 4.258 


m, breaker height \h, = 4.258 m, and breaker length = 
b 


0.692 dan —2 = 0.637. 


H 
Lp 


Breaking occurs when “= 0.. Therefore, from (14), 


breaking occurs when(1 = 5 = 0,, and subsequently, 


; ae ; kA 
there is a reduction in wave amplitude when ad 1. 


VI. REFRACTION - DIFFRACTION MODEL. 

In this section, a comprehensive two-dimensional wave 
transformation model is formulated, encompassing 
refraction, diffraction, shoaling, and breaking phenomena. 
The diffraction considered here is solely due to 
bathymetric features and does not account for the 
dispersion of wave energy in a direction perpendicular to 
the wave's propagation. 


The refraction-diffraction model was derived based on the 
shoaling-breaking equations elaborated in the preceding 
section. In the shoaling-breaking model, wave propagation 
occurs along the horizontal -x axis. However, in this 
section, the wave travels in the direction of the a, forming 
an angle a with the horizontal -x axis, as illustrated in 
Figure 2. 


Fig (2). Coordinate system, the relationship between č 
and x. 


The relationship between the € -axis and the x-axis is given 
by, 


E(x) = xcosa 


Where 
dé 
a TOSA arta (16) 
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For the wave number k which is a function of č, denoted 
as k = k(E(x)), 
dk dk dé 
dx dé dx 
substituting x in (13) with € and then replacing it in 7 ; 
while z is substituted with (16). 
dk 2k dh 
— = — — cosa 
dx 3A) dé 

(2n+>) 
Similarly, for the wave amplitude A = A(é(x)) 
dA dA 


dx dg" 
dA G (1 a) T 
dx E 5 cosh(@m) cosa 


where = is equation (16), replacing x to €. 


The calculation involves changes in wave direction. The 
wave direction a@ is 


v 
a = atan e) eee (17) 
da u2- 
bute dx dx 
i ape tti(‘é~*«C SSR (18) 


Where u represents the particle velocity in the horizontal- 
x direction, and v represents the particle velocity in the 
horizontal- y direction. To calculate these velocities, a 
three-dimensional velocity potential is utilized. 


(x,y) = Gcoshk(h + z) 


cos k(x cos æ + y sin æ) sin øt 


Where 

u(x, y,z,t) = a 

v(x, y,z,t) = -a 

4 The Results of Refraction-Diffraction Model. 


In the following section, an example of the execution 
results of the refraction-diffraction model is presented for 
a coastal bathymetry configuration in the form of a bay 
(Fig (3)), with waves having a wave period T = 7.15 
seconds and a wave amplitude in deep water Ag = 1.2 
meters. The incoming wave angle forms an angle Ay = 1.2 
. The calculation parameters used include the deep water 
coefficient 0 = 2.1 and y, = 1.4. 
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Ys. 


a=0° 


2004 


100 


/ 
“100 / 


200 


T I T j j T j | T 
0 100 200 300 400 500 600 700 800 900 


Fig (4). Wave height contours with an incoming wave 
angle a = 0° 


Fig (5). 3D wave height contours with an incoming wave 
angle a = 0° 


In Figures (4) and (5), the spread of wave energy towards 
the bay's sides is evident, indicating the model's effective 
simulation of wave refraction-diffraction. Furthermore, the 
model accurately simulates shoaling and breaking 
phenomena. 


Subsequently, the model results are showcased for the 
same wave with an incoming wave angle of a = 30° (Fig 
(6)). The outcomes are displayed in Fig (7) and Fig (8), 
demonstrating the model's proficient simulation of 
refraction, shoaling, and breaking processes. 
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Incoming wave 


Fig (6). Contour bathymetry with an incoming wave a=0°? 
angle a = 30° 


200 


100 


100 


T T T T T 
0 100 200 300 400 500 600 700 800 900 
500 


Fig (7). Wave height contours with an incoming wave 0d 72 X Ais 


ov 


angle a = 30° | ra 
3005 |) ///// 7 \ \ I 


200 4 


100 +4 


04 | : 
-100 4 | | 
-200 4 | | 
Fig (8). 3D wave height with an incoming wave angle 3007 \ 2 he ee A 
a = 30° Ñ E 2 
400 4 WSs Z + 
In the subsequent analysis, the model is applied to a 500 l m 
K . : i 0 100 200 300 400 500 
submerged island or sandbar (Fig (9)). Waves with a period Fig (11) Wave height contours 


of T = 6.2 seconds and an amplitude of Aọ = 0.9 m are 
utilized. The use of smaller waves is intentional; 
employing waves with larger amplitudes, such as Ay = 1.2 
meters, would lead to the island being submerged. 
Consequently, the wave height contours would flatten, 
devoid of any visible variation. 
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ee 


Fig (12) 3D wave height 


The model execution results on the submerged island (Fig 
(11) and Fig (12)) also demonstrate that the model can 
simulate the phenomena of refraction-diffraction, shoaling, 
and breaking effectively. 


VII. CONCLUSIONS 


By solving the Laplace equation along the sloping 
bottom, conservation equations governing water waves are 
derived, including the wave number conservation equation 
and the energy conservation equation. These equations 
regulate variations in wave constants as waves transition 
from deeper to shallower waters. 


Utilizing these conservation equations, formulating the 
equation for alterations in water wave constants as waves 
progress from deep to shallow waters becomes 
straightforward. 


Alterations in wave constants are interconnected, with 
changes in wave amplitude influencing wave length and 
vice versa. Consequently, modifications in these wave 
constants must be addressed simultaneously. 


The outcomes from the 1-D shoaling breaking model and 
2-D refraction-diffraction model simulations reveal the 
model's ability to effectively replicate shoaling-breaking 
and _ refraction-diffraction phenomena influenced by 
bathymetry. A necessary advancement involves 
incorporating the diffraction phenomenon into the 2-D 
model, accounting for the spread of wave energy 
perpendicular to the wave's direction. 
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